Mon. Not. R. Astron. Soc. 000, [T]|S] (2004) Printed 7 February 2008 (MN WF^K style file v2.2) 



Equilibrium Configurations of Homogeneous Fluids in 
General Relativity 

M. Ansorg*, T. Fischer, A. Kleinwachter, R. Meinel**, D. Petroff and K. Schobel 

Theoretisch-Physikalisches Institut, University of Jena, Max-Wien-Platz 1, 07743 Jena, Germany 
7 February 2008 

ABSTRACT 

By means of a highly accurate, multi-domain, pseudo-spectral method, we investigate 
the solution space of uniformly rotating, homogeneous and axisymmctric relativistic 
fluid bodies. It turns out that this space can be divided up into classes of solutions. 
In this paper, we present two new classes including relativistic core-ring and two- 
ring solutions. Combining our knowledge of the first four classes with post-Newtonian 
results and the Newtonian portion of the first ten classes, we present the qualitative 
behaviour of the entire relativistic solution space. The Newtonian disc limit can only 
be reached by going through infinitely many of the aforementioned classes. Only once 
this limiting process has been consummated, can one proceed again into the relativistic 
regime and arrive at the analytically known relativistic disc of dust. 
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1 INTRODUCTION 

Self-gravitating bodies of constant density have always 
played a central role in the physics of gravitation. Contri- 
butions that have been most significant to the aspects of 
this subject that will be relevant to this paper can be di- 
■ vide d into (1) work done within Newton's theory of gravita - 
tiomlMaclaurir] (118011) .lPoincareHl885tl .iDvsonl Il89lll893t) . 

iLichtensteinl <1933ft . lWond < 19741) and lEriguchi fc Sugimotd 
Jl98ll) and (2) work done within Einstein's the- 
ory o f gravitation: ISchw arzschil dl <ll9ld). IChandrasekharl 
Jl967l). iBardeenl Il97ll). iButterworth fc Ipserl Jl97rJ) and 

iGonde k-Roshi ska fc Gourgoulhonl j2002l) . Despite this great 
investment of effort, it has not yet been possible to ex- 
plore the whole spectrum of Newtonian, let alone relativis- 
tic solutions, even if one restricts oneself to axial symmetry 
and stationarity. Using sophisticated computer programs, 
we believe that a great step can be taken toward painting 
the complete relativistic picture of uniformly rotating, ho- 
mogeneous and axisymmetric bodies. However, because of 
the many intricate details in this picture, particularly in 
the vicinity of limiting configurations, it is necessary to use 
a computer program robust enough to be able to render 
such limiting configurations and accurate enough to be able 
to distinguish between neighbouring solutions of Einstein's 
equations. For our investigations of homogeneous fluids, we 
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used a code based on multi-d omain pseudo-spectral methods 
jAnsorg et alJ l200l l2003al) . which satisfies th ese require- 
ments . For a comparison with other codes, see IStergioulasI 

It is our intention in this paper to present the relativis- 
tic picture in its entirety by describing those parts of it that 
we have studied explicitly and augmenting them with con- 
jectures as to the rest. With this in mind, we focus our at- 
tention neither on the numerical methods nor primarily on 
properties of individual configu rations (as in lAnsorg et alJ 
l2003bl : fSchobel fc Ansordl2003l) . but instead emphasize the 
interrelations between various configurations and portray 
the solution space as a whole. As an interesting example 
of the newly explored configurations, we present the shape 
and various physical parameters of a core-ring solution. 

The picture that emerges for relativistic homogeneous 
fluids contains familiar demarcations, which we can use to 
orient ourselves. It contains, for example, three analytically 
know n solutions: the (inn er and outer) Schwarzschild solu- 
tion ( Schwarzschild 1916), the relativistic disc of dust solu- 
tion (iBardeen fc Wagonerll97ll:fNeugebauer fc Meinell99l 
2003) as well as the Maclaurin solution jMaclaurinl Il80 J) 
as one of its Newtonian limits. This last solution will be a 
particularly useful point of departure for describing the cor- 
responding relativistic picture. It represents homogeneous, 
rotating fluid spheroids and, but for a scaling factor, de- 
pends for given mass density on only one parameter, thereby 
allowing one to refer to the "Maclaurin sequence" . After in- 
troducing some basic equations in Sec. |2] we thus turn our 
attention to the Maclaurin solution and its post-Newtonian 
extension in Sec. [3] It turns out that not all relativistic con- 
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figurations are connected to each other in a continuous way, 
and it will be useful to introduce the notion of a class of 
solutions. A given class will be defined to include all config- 
urations of strictly positive mass that are connected contin- 
uously to each other. In Sec. [1] we review the characteristics 
of the "first" three classes in detail and provide an overview 
of the remaining classes. We close with a discussion of the 
limitations of numerical methods and some thoughts on the 
completeness of the relativistic picture painted in this paper. 



2 EINSTEIN'S FIELD EQUATIONS AND 
THEIR NUMERICAL SOLUTION 

Using Lewis-Papapetrou coordinates, the line element for a 
stationary, axially symmetric spacetime corresponding to a 
rigidly rotating perfect fluid configuration can be written in 
the form 

ds 2 = e 2a (dg 2 + d( 2 ) + W 2 e~ 2v {dip - to dt) 2 - e 2u dt 2 , 

where the metric functions a, W , v and uj are functions of g 
and £ alone. The boundary of a fluid configuration rotating 
with the angular velocity Q. obeys the equation 

e 2v - W 2 (lo - Q.) 2 e~ 2 " = const. = (f + Z y 2 , 

where Zq is the relative redshift observed at infinity of a 
photon emitted with zero angular momentum from the con- 
figuration's surface. Together with regularity requirements 
along the axis of rotation, g — 0, and the asymptotic be- 
haviour, the boundary equation including transition condi- 
tions and the field equations themsel ves form a complete 
set o f equations to be solved (see e.g. iButterworth fc Ipserl 
Il976l) . In this paper we choose natural units in which the 
gravitational constant G and the speed of light c are both 
equal to one. 

This set of equations was solved numerically to 
very high accuracy using a multi-domain pseudo-spectral 
method. In the algorithm used, the domains are chosen such 
that the unknown boundary £ = ( s (g) of the configura- 
tion always coincides with the boundary between two do- 
mains, thus avoiding Gibbs' phenomena. The boundary en- 
ters into the equations through coordinate transformations, 
which map the physical domains onto a square. The metric 
functions and the boundary are represented by a truncated 
Chebyshev expansion so that the above mentioned set of 
equations can be reduced to a set of (non-linear) algebraic 
equations. Typically, we used between f and 30 Chebyshev 
coefficients per dimension, depending on the properties of 
the function being represented. The parameters we can set 
to specify a fluid configuration include a mass-shedding pa- 
rameter (/3 of Sec. 14. in and central pressure, which can be set 
strictly to infinity (through a rescaling) , thereby allowing us 
to study even the boundary configurations to be discussed in 
th is paper with g reat ac curacy. Further details can be found 
in lAnsorg et alJ J2003al) . 

3 THE MACLAURIN SEQUENCE AND ITS 
BIFURCATION POINTS 

Maclaurin showed that for a spheroid (also called an ellipsoid 
of rotation) of constant density there is a balance of gravita- 
tional, centrifugal and pressure forces, leading to equilibrium 



Table 1 . Values for the ratio of polar to equatorial radius of the 
Maclaurin spheroid at the first six bifurcation points associated 
with the onset of axially symmetric, secular instability. 



Bifurcation points 



Ai 


= 0.17126186 . . . 


A 2 


= 0.11160323 . . . 


A 3 


= 0.08274993 . . . 


A 4 


= 0.06574427 . . . 


A 5 


= 0.05453402 . . . 


A 6 


= 0.04658868 . . . 



in Newtonian theory. For a given mass density, this spheroid 
can be described entirely by specifying the ratio of polar to 
equatorial radius, 

A := n/ra, (1) 

and the focal distance of the ellipse used to generate this 
ellipsoid of revolution. The focal distance merely provides 
the scaling of the spheroid and we can use A £ [0, 1] to 
par ametrize t he w hole s olution. 

iPoincarel (Il885l). iBardeenl dl97lD and 
iHa^m^u^^Eriguchil Jl984h all studied points of insta- 
bility along the Maclaurin sequence and the latter three 
authors were able to show that these correspond to bifur- 
cation points. The bifurcation points that are of particular 
interest to us here are those marking the onset of the 
various modes of axially symmetric, secular instability. 
Values for A at the first few of these can be found in 
Tableland will be denoted by A n , (n — 1 , 2, 3, . . .). There 
are countably infinitely many such bifurcation poin t s with 
the accumulation point A = 0. IChandrasekharl <ll967f) 
calculated the post-Newtonian corrections to the Maclaurin 
spheroids and found t hat they becom e singular at the 
bifurcation point A\. IBardee conjectured that 

the n th post-Newtonian approximation to the Maclaurin 
spheroids becomes singular at each of the first n such 
bifurca t ion po ints - a conjecture which was then proved in 
IPetrofj i2003l) . 

Assuming the convergence of the post-Newtonian series, 
an assumptio n based on strong evid ence especially in the 
disc limit fsee lPetroff fc MeineJl200lh . the singularities then 
represent impermeable barriers in an appropriate parameter 
space. These barriers allow us to divide our solution space 
into "classes" of configurations in the following way: We first 
divide up the Maclaurin sequence into an infinite number 
of segments with the above mentioned bifurcation points, 
An, separating one segment from the next. These segments 
can be numbered (or named) starting, for example, from 
A = 1. We can then identify uniquely any configuration of 
strictly positive mass with a given segement by requiring 
that the segment be approachable in a continuous way. All 
those configurations that have been thus associated with 
segment n make up class n. 

As we shall see in the following sections, numerical re- 
sults lend strong support to the idea that such impermeable 
barriers as were mentioned above exist and that the defi- 
nition of a class is thereby justified. In particular, we shall 
see that it is possible to place ourselves at one of the end- 
points of a Maclaurin segment and move off the Maclaurin 
sequence along other well-defined boundaries until we land 
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again on the Maclaurin segment at the second end-point. 
Neighbouring classes have only this one point as a common 
boundary. 



4 THE RELATIVISTIC CLASSES 

4.1 Class I: The Generalized Schwarzschild Class 

The first class of configurations t hat will be discussed in thi s 
paper was explored extensively in lSchobel fc Ansord {2003). 
Here we shall give a brief overview of it, paying particular 
attention to its defining boundaries. 

For a given constant mass density, /i, a stationary, ax- 
isymmetric, rigidly rotating relativistic fluid depends on two 
parameters. A limiting configuration, i.e. one that possesses 
a physical property representing a boundary, is specified by 
fixing a corresponding parameter, meaning that it depends 
on only one further parameter. It is therefore possible to 
refer to a boundary sequence and identify the end-points of 
this sequence uniquely with respect to the second parameter. 

Clearly the static boundary is marked by the specifica- 
tion A = 1. (Note that we also use the radius ratio as defined 
in Eq. Q in the relativistic context, where ri and r-z are now 
coordinate radii.) Static configurations are described by the 
analytic Schwarzschild solution, which can be parametrized 
by their central pressure, p c . Starting from the static end of 
the Maclaurin sequence (A — l,p c — 0), let us increase the 
central pressure until it becomes infinite. This sequence cor- 
responds to a mere point in the lower right corner of Fig. Q 
and to the segment from (a) to (b) in the schematic diagram 
of Fig. 2] Configurations with infinite central pressure are 
also clearly boundary configurations. We can parametrize 
these by the normalized angular velocity and follow the se- 
quence from the static limit (where the angular velocity is 
of course zero) up to its maximal value marked by a mass- 
shedding limit ((b) to (c) in Fig. H^. Configurations at the 
mass-shedding limit rotate with an angular velocity equal 
to that of a test mass rotating at the equator on a circular 
orbit. They exhibit a cusp at the equatorial rim and are thus 
easy to identify by the geometry o f their surface . We de fine 
a mass-shedding parameter j3 as in lAnsorg et alJ <l2003al) by 



P 



— llm C< 



dg ' 



where £ = C^io) represents the boundary of the configura- 
tion (see Section The mass-shedding parameter takes on 
the value (3 = in the mass-shedding limit and j3 = 1 for 
Maclaurin spheroids. Fixing the value j3 — 0, we proceed 
along the mass-shedding sequence by reducing the central 
pressure until we reach the Newtonian limit p c — ((c) to 
(d))- 

We can trace the last part of this sequence in the en- 
largement in Fig. i.e. follow the red curve as it comes in 
from the top. This curve crosses the Maclaurin sequence in 
the parameters plotted there at a value for the radius ra- 
tio of A ~ 0.195. The two configurations at this cross-over 
point are however clearly not the same since one is a rel- 
ativistic star rotating at the mass-shedding limit whereas 
the other belongs to the Newtonian Maclaurin sequence. 
In some other parameter diagram these two configurations 
would correspond to distinct points. This and various other 
intersections in Fig. are merely a consequence of the fact 



that the configurations considered in this paper depend on 
two parameters, but not always in a unique way. 

We resume our journey along the boundary curves by 
following the Newtonian sequence, Af , from (d) to (e) 
(parametrized again by p c — for example) by increasing (3 
until we reach the Maclaurin curve at the value f3 = 1 and 
find ourselves precisely at the bifurcation point A = Ai. This 
Newtonian sequence and indeed all Newtonian sequences bi- 
furcating from the Maclaurin sequen ce at the first ten A n 
were studied in lAnsorg et all J2003cft in which a depiction 
of the shape of the configurations along these sequences was 
also provided. 

Except for the Schwarzschild solution itself, we consider 
all the boundaries in the classes to be open, i.e. the limit- 
ing sequences themselves are not considered to belong to 
a solution class. In the case of the Newtonian limit or the 
limit of configurations with infinite central pressure, it seems 
appropriate to exclude such configurations from the physi- 
cally permissable solutions. In the case of the mass-shedding 
limit or the two-body limit that will arise in later classes, 
this choice is a mere matter of convention. The extreme Kerr 
Black Hole, which will present itself in the next section, is 
also taken to be an open boundary. 

Having circumscribed the generalized Schwarzschild 
class by specifying its boundary configurations, it is pos- 
sible to discuss those configurations making up this class. 
Amongst the most interesting properties of this first class is 
that there is an upper limit on the attainable mass, which is 
much higher than that known from the static solution. The 
maximal gravitational mass M m 0.19435/^//! is roughly 
34% greater than for the corresponding non-rotating config- 
uration and is reached at (c). 

4.2 Class II: The Generalized Dyson Ring Class 

A portion of the gene ralized Dyson ring class was explored in 
lAnsorg et all (l2003bD . Here we extend that work and define 
the entire class by navigating its boundaries as with class I. 
We begin by placing ourselves on the Maclaurin curve at 
the point A = A\ and proceeding along the sequence A~[ by 
decreasing the value of A for My/Jl = 0. It turns out that 
configurations along this sequence pinch together in the cen- 
tre (A = 0), marking the transition from a sphero idal to a 
toroid al topology. This behaviour was predicted bv lBardeenl 
Jl9TlT) . who also correctly surmised that ring configurations 
are connected to the Mac laurin seq uence via the bifurca- 
tion p oint Ai . IWoiiel il974h as well as lEriguchi fc Sugimotol 
(1981) studied such configurations and their connection to 
the Maclaurin seque nce, aspects of which were clarified in 
lAnsorg et all j2003ch . 

If we define A for toroidal topologies to be the negative 
ratio of the inner to outer coordinate radius of the ring, i.e. 



A:=H< 

T2 



n 



'■■z 



polar radius 
equatorial radius 

-ri = inner radius 
= outer radius 



spheroidal 



toroidal 



then we can continue to decrease its value until we finally 
arrive at A — — 1. This limit, abou t which Dyson expanded 
his ring solution l|Dvsorll 892. 1893), is problematic however, 
since the "corotating Newtonian surface potential" tends to 
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A 




Figure 2. A three-dimensional parameter diagram depicting all of class III. Lines of constant M^fjl beginning with the Newtonian limit 
M^ffl = (green and blue curves) are drawn in steps of 0.01. The colour coding of the boundary sequences corresponds to that of Fig. HI 
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minus infinity if the Newtonian mass and the radius T2 are 
chosen to be strictly positive and finite. The open circle at 
(f) in Figs Q and 0] indicates the subtle and interesting jump 
to the extreme Kerr Black Hole limit (Zq — > oo). 

Although the horizon of the extreme Kerr Black Hole 
viewed from the "outside" world is a mere point at the 
origin of the coordinate system used here, there exists an 
"in ner" world (if the limit is taken in rescaled coordinates, 
see lAnsorg et aP feooab 5 iMeinell 1200211 in which we can in- 
crease the value of A until we come to a mass-shedding limit 
((f) to (g) in Fig.0. Following the mass-shedding sequence 
by decreasing the relative redshift, Zo, we cross back over 
to spheroidal topologies and arrive at the Newtonian mass- 
shedding configuration for a redshift of zero ((g) to (h)). Fol- 
lowing the Newtonian sequence by increasing the value 
of the mass-shedding parameter, /3, we arrive at the bifurca- 
tion point A2 of the Maclaurin sequence and have completely 
circumscribed the generalized Dyson- Ring class ((h) to (i)). 

The configurations in this class, in contrast to those of 
the generalized Schwarzschild class, contain no upper bound 
for the mass. This surprising fact would only be of astrophys- 
ical relevance if ring configurations turned out to be stable. 1 

Another salient feature of this class is that a continuous 
transition to an (extreme Kerr) Black Hole is possible. This 
can be compared to the first class for which configurations 
reached the boundary of infinite central pressure for finite 
values of the redshift, which has a global maximum of Zq « 
7.37825 in the generalized Schwarzschild class. Note that the 
extreme Kerr Black Hole is the only candidate for a Black 
Hole limit of rotating fluid bodies in equilibrium llMeinell 
|2004) . 

4.3 Class III: The First Generalized Core-Ring 
Class 

In this section, we turn our attention to the hitherto un- 
explored class III and locate its boundaries, depict it in a 
three-dimensional parameter space and provide a represen- 
tative example of one of its configurations. Having arrived 
at A = A2 from class II by increasing the value of A (cf. 
Fig. , we proceed into the third class along the Newtonian 
sequence A\~ by placing ourselves at A = A2 and increas- 
ing the value of this ratio further. The configurations take 
on a central bulge shape, but begin to pinch in toward the 
outer edge. Finally the indentation pinches together yield- 
ing a ring just barely in contact with a central core corre- 
sponding to the pat h from (i) to (j) in Fig. 0] (see Fig. 7 in 
lAnsorg et aTll2003cJ) . Although it would be possible to con- 
sider the separation of core and ring, we are not interested 
here in studying a many-body problem and restrict our at- 
tention to single, homogeneous bodies. Thus this "pinching 
together" marks an end-point of the Newtonian sequence. 
We can follow the sequence of configurations on the verge of 
forming a two-body system by allowing the mass to increase. 
This sequence ends in a mass-shedding limit for the configu- 
ration with a global maximum for M^JJl m 0.16 in this class 
((j) to (k) in Fig.^J. Placing ourselves on the mass-shedding 
sequence, we can allow the mass to decrease until we reach 

We intend to investigate the stability of such configurations in 
the future, but preliminarily suppose that they are unstable. 



Table 2. Values for dimensionless physical quantities for an ex- 
emplary configuration in Class III with increasing order, m (see 
lAnsorg et all2003al) . of the numerical approximation. The quan- 
tities not yet referred to in this paper, Mo, J and -R c i rc are bary- 
onic mass, angular momentum and circumferential radius respec- 
tively. The position of this configuration in Figs0and|5]is marked 
by the letter a. 





m = 20 


m = 30 


m = 40 


A 


0.1 


0.1 


0.1 




0.1 


0.1 


0.1 


M 0v /7i 


0.1093062 


0.10930610 


0.109306099 


Jfi 


0.01483899 


0.014838966 


0.0148389661 




0.8627113 


0.86271067 


0.862710640 


^circV^ 


0.5843826 


0.58437986 


0.584379882 




0.4921657 


0.49216554 


0.492165525 



C/T2 




Q/r2 



Figure 3. The meridional cross-section of the boundary of the 
Class III configuration described in Table l2l 

the Newtonian limit ((k) to (1)). From here, the Newtonian 
sequence can be followed to the point A = A3 by increasing 
the parameter j3 ((1) to (m)). 

As was discussed in Sec. 14.11 two parameters do not 
necessarily suffice to describe a configuration uniquely. By 
adding a third dimension, one can resolve such ambiguities. 
In Fig. [5] a plot of the A, f3 and M^fji values for the config- 
urations of Class III can be found. In this plot, a given point 
on the two-dimensional surface in the three-dimensional pa- 
rameter space refers only to a single configuration. A projec- 
tion of the boundary curves onto the ^4-/3 plane would yield 
a picture of similar complexity to the depiction of class III 
in Fig. [I] 

An example of various physical parameters for a con- 
figuration with the prescribed values A = 0.1, M-JJi — 0.1 
can be found in Table [5] The shape of the boundary of this 
configuration in meridional cross-section is shown in Fig. 

4.4 The Further Classes 

As one moves to higher and higher classes, the configurations 
tend to grow more disc-like. The distance from the Maclau- 
rin curve in the parameters of Fig. Q grows small although 
the shapes of the configurations deviate significantly from a 
spheroid. We term the odd classes, beginning with the third, 
"core-ring" classes because the configurations on the verge 
of forming a two-body system contain a central body with 
a surrounding ring. The central body will itself grow flat in 
shape and corrugated in appearance as one moves to higher 
classes. Beginning with the fourth class, the even ones, on 
the other hand, tend to grow concave at the centre as they 
"emerge" from the bifurcation points. They then pinch to- 
gether at the centre and proceed into the toroidal region in 
which A is negative. Here too, the ring configurations grow 
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flat and corrugated as one moves to higher classes. The cor- 
rugations become more pronounced until finally one of them 
(presumably the outermost) pinches together, whence the 
name "two-ring" classes for such configurations. The num- 
be r of corrugations of t he Newtonian sequences is described 
in lAnsorg et al.l l)2003cT) . There it was shown at least up to 
A 10 that the two Newtonian sequences bifurcating from each 
A n end in a mass-shedding limit and a two-body configura- 
tion respectively. We assume that this holds for all the bi- 
furcation points. Emanating from the two bifurcation points 
associated with a given class are two Newtonian sequences, 
the one ending in a mass-shedding limit and the other in a 
two-body system. Insofar as the Black Hole limit and that 
of infinite central pressure are precluded from the higher 
classes (as indicated by our numerical investigations), the 
two end-points of the Newtonian sequences are connected 
by a mass-shedding and a two-body sequence. The point at 
which these two configurations meet most likely marks the 
configuration with the greatest mass for that class. 

In order to arrive at this picture for the higher classes, 
we made use of the first ten Newto nian bifurcation sequences 
presented in lAnsorg et aD d2003cl) and the post-Newtonian 
behaviour in order to make a conjecture, which could then 
be verified by studying two new classes (III and IV). The 
properties of the remaining classes could then be extrapo- 
lated and, as such, cannot be considered absolutely certain, 
but rest on firm ground. 

In Fig. [1] the classes are depicted schematically. As in 
this depiction, two neighbouring classes have only the bifur- 
cation point along the Maclaurin sequence in common. The 
first two classes are atypical in that the first contains the 
static boundary sequence and one with infinite central pres- 
sure and the second contains a sequence of infinite redshift 
(extreme Kerr Black Hole sequence). Furthermore, these are 
the only two classes in which we found ergo-regions. Our 
investigations show that they are always toroidal in topol- 
ogy and "begin" to appear within the fluid configuration. A 
more detailed discuss ion of the ergo-regions can be found in 
I Ansorg et aD (l2003allbh . 

The definition of the classes given in Section [3] relies 
on the fact that neighbouring classes have only the bifurca- 
tion point along the Maclaurin sequence A n as a common 
boundary. As was already mentioned, the post-Newtonian 
approximation strongly suggests that any sequence of con- 
stant, strictly positive mass, parametrized by /3 and A for 
example, will be deflected away from the singular point 
ri/r2 = A n and "back into" the original class. This is ex- 
actly what is observed numerically. In Fig. "5J /3 versus A 
values are plotted for sequences of constant mass. Numer- 
ical investigations show that sequences belonging to differ- 
ent classes are completely disjoint although the gap in the 
vicinity of a given A n grows small for decreasing mass and 
approaches the poin t (A n ,l)- This feature, which was al- 
ready conjectured bv lBarde'enl (Il97lf) . can be seen distinctly 
in the figure around A±. Such a gap was observed numeri- 
cally around A2 and A3 as well. A comparison with Fig. [2] 
shows why the curve M^/JI = 0.01 of that plot deviates so 
markedly from the Newtonian boundary curves. If one were 
to plot such contour lines for increasingly diminishing val- 
ues of My/JI, the tendency to approach the Newtonian limit 
would be evident. 




0.14 0.16 Ai 0.18 0.20 

A 

Figure 5. The relationship of /3 to A for sequences of constant 
mass in the vicinity of the bifurcation point associated with A\ 
(the point (e), cf. Fig. |1J. The solid line represents configurations 
with My/jl = 10 -4 , the dashed line with M^JJl = 10~ 5 and the 
dotted line with Myfji = 10 -6 . The thin horizontal line represents 
the Maclaurin sequence and the thin line crossing it is made up 
of the Newtonian bifurcation sequences and A~^ (see Fig, 14*1. 
Note the gap between the sequences of class I and class II. 

5 DISCUSSION 

Numerical evidence suggests that as one proceeds to flat- 
ter configurations (higher classes), the departure from the 
Newtonian sequence grows small. This is reflected in the 
fact that boundaries associated with highly relativistic at- 
tributes such as infinite central pressure or the transition to 
a Black Hole are presumably absent from all classes upwards 
of class II and is a result of excluding many-body configu- 
rations. If one were to abandon this restriction and choose 
some (arbitrary) relationship for the rotation of the various 
segments, then it would again be possible to reach "relativis- 
tic boundaries" most likely. With this restriction in place, 
however, the only path to the disc of dust passes through an 
infinite number of bifurcation points and sees the configu- 
rations growing ever flatter and more and more corrugated 
until one lands necessarily at the Newtonian Maclaurin disc. 
Only once this limit has been reached can one proceed again 
into the relativistic regime, i.e. the relativistic disc of dust 
iNeugebauer fc Meinellll995li and indeed reach the extreme 
Kerr Black Hole. Th is picture reinforces Bardeen's specula- 
tions jBardeerJl97lTl . who wrote that the singularities in the 
post-Newtonian expansion at the bifurcation points "may 
forbid the existence of any highly relativistic, highly flat- 
tened, uniformly rotating configurations which are simply 
connected." He went on to conclude that the "only accept- 
able model of an infinitesimally thin, uniformly rotating disk 
in general relativity may be one that is made up of an infinite 
number of disjoint rings." 

The complicated and highly non-linear set of equations 
describing axially symmetric, stationary, uniformly rotating 
fluids of constant density can only be solved approximately 
or numerically. As such, it is difficult to find strict results 
describing the solution set. It is conceivable, for example, 
that there exist solutions to these equations that possess no 
Newtonian limit (and must thus be disconnected from the 
classes of solutions introduced here). It is also conceivable 
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Figure 4. A schematic portrayal of the classes and their boundaries. A depiction of the surface of configurations in meridional cross- 
section is provided for configurations at the intersection point of two boundary sequences (except for the point (f), see Sec. 14.21 . In 
addition, various surfaces are shown to help illustrate the transition from one corner to the next. Whenever ergo-regions arise, they are 
depicted by dashed lines (hence the dot in the centre of the configuration located at (b)). 
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that Newtonian solutions exist that have not yet been found, 
but serve as limits to relativistic solutions. 

Although our numerical considerations do not preclude 
the possibility of the existence of undiscovered solutions, we 
presume to conjecture that the classes presented here cover 
the entire solution space. Attempts to prove or disprove this 
conjecture are bound to lead to innovative insights into the 
nature of Einstein's equations and novel methods for probing 
their structure. The modification to the picture drawn in this 
paper resulting from a change in the equation of state will 
be discussed elsewhere. 
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